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Context. Stellar jets are an essential ingredient of the star formation process and a wealth of information can be derived from their 
(-H ' characteristic emission-line spectra. 

Aims. We investigate the physical structure and dust reprocessing in the shocks along the beam of a number of classical Herbig-Haro 
(HH) jets in the Orion and Lupus molecular cloud (HH 1 1 1, HH 1/2, HH 83, HH 24 M/A/E/C, Sz68). The parameters describing the 
plasma conditions, as well as the dust content, are derived as a function of the distance from the source and, for HH 1 1 1, of the gas 
velocity. 

Methods. Spectral diagnostic techniques are applied to obtain the jet physical conditions (the electron and total density, n e and n H , 
the ionisation fraction, x c , and the temperature, T f ) from the ratios between selected forbidden lines. The presence of dust grains is 
investigated by estimating the gas-phase abundance of calcium with respect to its Solar value. 

Results. We find the electron density varies between 0.05-4 10 3 cm" 3 , the ionisation fraction x f is 0.01-0.7, the temperature ranges 
between 0.6-3 10 4 K, and the hydrogen density between 0.01-6 10 4 cnT 3 . Interestingly, in the HH 111 jet, n e , x e , and T c peak in 
the High Velocity Interval (HVI) of the strongest working surfaces, confirming the prediction from shocks models. Calcium turns out 
to be depleted with respect to its Solar value, but its gas-phase abundance is higher than that estimated in the interstellar medium in 
Orion. The depletion is high (up to 80%) along the low-excited jets, while low or no depletion is measured in the jets which show 
higher excitation conditions. Moreover, in HH 111 the depletion is lower in the HVI of the faster shock. 

Conclusions. Our results confirm the shock structure predicted by models and indicate that the shocks occurring along the jets, and 
presumably those present in the launch zone, are partially destroying the dust grains and that the efficiency of dust reprocessing 
strongly depend on shock velocity. However, the high Ca gas-phase abundance estimated in some of the knots, is not well justified by 
^\ t existing models of dust reprocessing in shocks, and indicates that the dust must have been partially reprocessed in the region where 

f*"^ . the flow originates. 
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' 1. Introduction lighting the stratification in the shock cooling regions (Nisini 

et al. 2005; Podio et al. 2006), the variation of the excitation 

Although it is widely accepted that jets play an essential role in conditions along the jet (Bacciotti & Eisloffel 1999), across it 

star formation, there are still many open questions about how (Bacciotti et al. 2000; Hartigan & Morse 2007) and, in a few 

they are generated and propagate in the interstellar medium, cases, with gas velocity (Lavalley-Fouquet et al. 2000; Coffey et 

Their optical spectra are characterised by emission lines from al. 2008; Garcia Lopez et al. 2008). 

atomic species which are collisionally excited in the shock An important aspect is also the estimate of the dust content 

waves generated by the interaction of the ejected material with in jets, a topic that has been poorly investigated to date. Several 

the interstellar medium or previously ejected jet material. These theoretical studies have explored the dust content and disnibu- 

lines contain a mine of information on the jet and the shock prop- tion in the interstellar medium. On one hand both observations 

erties. Beyond the jet/shock morphology and kinematics, infor- and theoretical models show that the gas-phase abundances of 

mation on the gas physical conditions can be derived by devel- elements like iron (Fe), magnesium (Mg), silicon (Si), and cal- 

oping a grid of shock models the predictions of which are com- cium (Ca) are considerably depleted (by a factor of the order 

pared with the observed line ratios (e.g., Raga & Bohm 1986; of 10 2 -10 4 ) with respect to their Solar abundances because of 

Hartigan et al. 1987, 1994); or, alternatively, by using spectral di- the formation of dust grains (Baldwin et al. 1991; van Dishoek 

agnostics techniques which are independent of any assumption et al. 1993; Savage & Sembach 1996). On the other hand the 

on the heating mechanism (e.g., Bacciotti & Eisloffel 1999). The chemical composition, the size and the content of dust grains 

application of these diagnostic techniques to various datasets al- can evolve as they lose atoms to the gas phase because of shocks 

lowed derivation of the jet physical/dynamical structure, high- waves. The high energy gas-grain and grain-grain collisions oc- 
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curring behind the shock front, in fact, may lead to the erosion 
of the grain surfaces (sputtering) and/or to the vaporisation and 
fragmentation of grains (Jones et al. 1994; Jones 2000; Draine 
2003; Guillet et al. 2009, and references therein). Total dust de- 
struction is expected in high velocity shocks such as supernova- 
generated shock waves (v s > 200 km s" 1 ) but there are very few 
observations about dust reprocessing in slower shocks such as 
the ones occurring along stellar jets (v s ~10-80 kms~'in the jet 
beam). A few works have investigated the gas phase abundance 
of Fe in HH jets (Beck-Winchatz et al. 1994, 1996; Mouri & 
Taniguchi 2000; Bohm & Matt 2001; Nisini et al. 2002) and 
only recently, the analysis has been extended to other refractory 
species such as Si, Ca, carbon (C), nickel (Ni), chromium (Cr), 
and titanium (Ti) (Nisini et al. 2005, 2007; Podio et al. 2006; 
Garcia Lopez et al. 2008). 

To further investigate these aspects, in this paper we study 
the distribution of physical parameters and, in particular, the gas 
phase abundance of calcium in a sample of HH jets (HH 111, 
HH 1/2, HH 83, HH 24 M/A/C/E, Sz68), as a function of the 
distance from the source, i.e. in the different working surfaces 
along the jet beam, and, where possible, of the gas velocity. As 
we will show, these kind of estimates are very important to un- 
derstand dust reprocessing in shocks and are also useful to derive 
constraints on the properties of the jet launching region. 

This paper is organized as follows: in Sect. 2 we present 
the observations and the data reduction process and in Sect. 3 
we briefly recall the principle of the so-called BE diagnostic 
technique, which has been used to derive the gas physical 
conditions along the jets. In Sect. 4 we present the physical 
properties derived for HH 1 1 1, in two velocity intervals, and for 
the jet from Sz68, that is analysed here for the first time. Similar 
results derived for the other jets of the sample at higher spatial 
sampling than in previous studies are presented in the online 
material. In Sect. 5 we use the obtained physical parameters to 
estimate the presence of dust grains in the jets and we discuss 
the efficiency of shocks in reprocessing the dust comparing our 
results with the predictions of theroretical models. Finally, in 
Sect. 6 we present a summary of our results. 



in the other lines may be off the slit center inducing a velocity 
shift in these lines (see appendix from Bacciotti et al. 2002 for a 
discussion of the uneven slit illumination effect). The resampling 
in velocity and the correction for the above mentioned velocity 
shifts allowed us to accurately align the lines in velocity despite 
the low spectral resolution. 

As shown in Hartigan & Morse (2007), the comparison be- 
tween the slit width and the projected width of the target is cru- 
cial in spectroscopic observations. If the slit is narrower than 
the jet (and the lines are broader than the spectral resolution) 
the detected lines are spectrally resolved and the velocity res- 
olution is higher for decreasing slit width. If, on the contrary, 
the slit is much wider than the jet (and the velocity dispersion 
is lower than the spectral resolution) the observations produce 
emission-line images for each line in the spectrum, but no ve- 
locity information can be recovered. The latter is the so-called 
"slitless spectroscopy" (see also Hartigan et al. 2004). 

Our observations are seeing-limited and the jet widths are 
spatially unresolved. Thus the critical parameter to be compared 
with the slit width is the Full Width Half Maximum (FWHM) 
of the seeing. This was varying between when observing 
HH 1/2, HH 83, HH 24 M/A/C/E, and Sz68, and ~1'.'5 when 
observing HH 111. 

In the case of HH 111, the seeing was comparable to the 
slit width and the acquired spectra show line profiles which are 
broader than the spectral resolution. Thus the HH 1 1 1 jet is par- 
tially resolved in velocity, allowing the analysis of its structure 
along the jet and in two different velocity intervals, as we will 
explain in more detail in the next section. 

In the other jets, on the contrary, the slit is wider than the 
seeing FWHM. Thus we obtain images of the jet in the differ- 
ent lines. This is evident when comparing spectral images of the 
detected lines with pre-imaging acquired with the [S n] filter. To 
avoid the confusion between the spectral and spatial informa- 
tion we integrated the line fluxes over the line "spectral" profile, 
suppressing the velocity dimension. 

In the next section we will explain in more detail the analysis 
performed in the two different cases. 



2. Observations and Data Reduction 

We observed a sample of HH jets (HH 111, HH 1/2, HH 83, 
HH 24 M/A/C/E, Sz68) in the optical range. Our spectra were 
acquired in January 1998 at the ESO 3.6m-telescope by using 
the EFOSC2 spectrograph, the GR9 grism and a 1'.'5 slit width. 
The slit was aligned parallel to the jet axis with position an- 
gles: 276.7° for HH 111, 324.7° for HH 1/2, 297.7° for HH 
83, 328.1° for HH 24 C/E, 219.0° for HH 24 G, and 334.9° for 
Sz68. The spectra cover the wavelength range from 5600 to 7335 
A with a spectral resolution -1300 (FWHM iM „.~230 kms -1 ). 
The spatial and spectral scales are 0. 157"/pixel and 0.872A/pixel 
corresponding to ~40 km s _1 . 

The spectral images were flat-fielded, sky-subtracted, wave- 
length and flux calibrated. In addition to the standard 
data reduction process the spectral images of the detected 
lines ([Oi]X16300,6363, [Nn]Xl6548,6583, [Sn]/Ll67 16,6731, 
[Can]/L1729 1,7324) were velocity-calibrated, resampled to the 
same velocity scale for all the lines, and corrected for velocity 
shifts due to atmospheric differential refraction. During the ob- 
servations, in fact, the slit was aligned along the jet emission 
through pre-imaging in the [Sn]/16731 line. Because of atmo- 
spheric differential refraction, however, the peak of the emission 



3. Application of Spectral Diagnostic 

Our analysis relies upon the dependence of the ratios between 
optical forbidden lines on the gas physical conditions in the jet, 
i.e. on the electron and total hydrogen density, n e and n#, the 
hydrogen ionisation fraction, x e , and the electron temperature, 
T e (Bacciotti et al. 1995; Bacciotti & Eisloffel 1999). 

As a first step, we thus computed for all the jets 
in our sample the line ratios used in the diagnostics, 
i.e. [Sn]/I6731/6716, [Nii]/L16548,6583/[Oi],U6300,6363, and 
[Oi]/i/i6300,6363/[Sii]/i/i6716,6731.ForHH 111, which is par- 
tially spectrally resolved, we divided the spectral images of the 
detected S + , O , and N + lines pixel by pixel, thus obtaining 
position-velocity (PV) diagrams of the line ratios (see Fig. 1). 
For the other jets (HH 1/2, HH 83, HH 24 M/A/C/E, Sz68) we 
integrated over the spectral profile obtaining the variation of the 
ratios only as a function of the distance from the source. The dis- 
tribution of the ratios already gives a qualitative idea of the vari- 
ations of the gas physical conditions along the jet and, for HH 
111, with the gas velocity. The ratio between the sulphur doublet 
lines ([S n]/1673 1/6716) is a tracer of the gas electron density up 
to the critical density of the [Sn] lines (n c ~2.5 10 4 cm -3 ), the 
[N n]/[0 1] increases primarily with x e , and the [O i]/[S n] ratio 
is dependent on both x e and T e . 
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To obtain a more quantitative estimate of the gas physi- 
cal properties we use the so-called BE technique (Bacciotti & 
Eisloffel 1999), that extracts diagnostic information from the 
comparison of observed and calculated ratios. The method is 
based on the fact that the involved forbidden lines are colli- 
sionally excited and that, in low-excitation conditions and pro- 
vided that no strong sources of ionizing photons are present, the 
ionization fraction of hydrogen is tightly related to those of ni- 
trogen and oxygen via charge exchange processes. This allows 
one to easily retrieve x e (and the total hydrogen density) and T e 
from the comparison of the observed line ratios. As discussed in 
Bacciotti & Eisloffel (1999) this method assumes that the lines 
used in the technique are emitted in a region of similar tem- 
perature, density and ionization fraction. Bacciotti & Eisloffel 
(1999) and Nisini et al. (2005) carefully verified this assump- 
tion by using the results from shock models (e.g., Hartigan et 
al. 1987). They showed that the considered forbidden lines have 
their peak emission in the same region behind shock fronts of 
different shock velocity, and verified that the obtained parame- 
ters are representative of the average value that those quantities 
assume behind the shock front in the line emission region. 

In this paper we use an improved version of the BE code in 
which the values of the collision strength are updated using the 
results in Keenan et al. (1996) for S + , Hudson & Bell (2005) 
for N + , and Berrington & Burke (1981) and Mendoza (1983) for 
O . This allows us to obtain a better fit of the collision strength 
for high temperatures (T e between 2-5 10 4 K) and thus to trace 
correctly the plasma physical conditions in regions of high ex- 
citation, i.e. where the [Nn] emission is comparable to that in 
the [Sn] and [Oi] lines, like in some knots of the HH 83 and 
the HH 24 jets where we could not obtain results in our previous 
work (Podio et al. 2006). 

Note that the diagnostic technique uses ratios between dif- 
ferent species and, thus, the obtained \ e and T e may depend on 
the chosen set of elemental abundances. However, as stressed 
in Podio et al. (2006), the relative variations of the parameters 
(i.e. in different knots and/or velocity intervals) do not depend 
on the choice of an abundance set (see Fig. 1 of Podio et al. 
2006) which, therefore, can be thought as a model parameter. 
On the other hand the most recent determinations by Asplund et 
al. (2005) (Solar) and Esteban et al. (2004) (Orion) are in good 
agreement and give rise to diagnostic results that differ of at most 
15%. One of the main goal of this paper is the determination of 
the calcium gas-phase abundance with respect to Solar. Thus, at 
difference with Podio et al. (2006), the Solar abundances esti- 
mated by Asplund et al. (2005) are adopted for consistency in 
all the diagnostic analysis. In particular, the values of the tem- 
perature and the ionization fraction are larger when using the 
abundances from Esteban et al. (2004). The uncertainty over the 
absolute values is not shown in Fig. 3, 4, A.l, A. 3, A. 4, because 
these are intended to be used to analyse the relative variations 
of the parameters along the jet and in different velocity inter- 
vals. However, we do consider the uncertainty due to abundances 
choice when computing the theoretical [Ca n]/[S n] ratios to es- 
timate the calcium gas-phase abundance (see Sect. 5 for more 
details). 

Since the observations are seeing-limited, i.e. the angular 
resolution of the data is much lower than the spatial sam- 
pling, we integrated the line fluxes over the seeing FWHM 
(FWHM S( , e! „£~1.2"~8 pixels) before applying the diagnostics. 
Thus the de facto spatial sampling of the obtained results in all 
the figures is equal to the seeing. Moreover, we integrated the 
line fluxes over the spectral profile for the spectrally unresolved 
jets (HH 1/2, HH 83, HH 24 M/A/C/E, Sz68) and over two veloc- 



ity intervals for HH 111: Low (for v>- 1 00 km s~ 1 ) and High (for 
v<-100 kms -1 ) Velocity Interval (LVI and HVI, respectively). 

For the HH 111 and the HH 1 jet we corrected the observed 
fluxes for reddening using the values of the visual extinction, Ay, 
estimated by Nisini et al. (2005) and Podio et al. (2006), the stan- 
dard dereddening procedure from Draine (1989) and an interpo- 
lation of the extinction law derived by Rieke & Lebofsky (1985) 
for the near-IR bands. For the HH 83 and the HH 24 jet we have 
no estimates of the visual extinction because the [Fen]/I1.64, 
1.32, 1.25 lines were not detected in our previous optical/NIR 
spectra (Podio et al. 2006). Thus we could not correct the line 
fluxes for the reddening before to compute the ratios and apply 
the diagnostics. On the other hand, Bacciotti & Eisloffel (1999) 
showed that, since the lines used in this diagnostic are very close 
in wavelength, the errors induced by the fact that we are not cor- 
recting for extinction are at most about 8-10% for the ionization 
fraction and about 15% for the temperature as long as the visual 
extinction is lower than ~3. 

Once a set of abundances is assumed, the errors affecting 
the parameters obtained through the BE technique are due to the 
signal-to-noise ratio of the measured fluxes and the uncertainty 
in the determination of the Ay which is used to deredden the line 
fluxes. These are very low (<5%) in the bright knots of the HH 
111, HH 1, and HH 83 jets which have strong emission in all the 
lines, and larger (up to 50%) in the fainter knots of the HH 24 
and Sz68 jets. 



4. The Physical Structure of the Jets 

In this section we briefly describe the new information gathered 
by applying spectral diagnostics (i.e. the BE technique) to our 
3.6m/EFOSC2 spectra of a sample of HH jets. 

In Sect. 4.1 we present the physical parameters (n e , x e , T e , 
Tin) for HH 111, for which we have been able to highlight for 
the first time the stratification of the gas physical conditions as a 
function of velocity. 

In Sect. 4.2 we present the parameters obtained for the HH 
jet from Sz68 which has not been analysed before with the BE 
technique. 

The other objects in our sample, already analysed previously 
(Bacciotti & Eisloffel 1999; Nisini et al. 2005; Podio et al. 2006), 
are here re-analysed with the improved version of the code and 
increased spatial sampling, for an accurate determination of Ca 
gas-phase abundance along the jets (see Sect. 5). The physical 
properties of these jets are shown in the online material, for ref- 
erence. 

4.1. HH 111: the High and Low Velocity Intervals 

HH 111 is a well-known parsec-scale jet powered by the young 
star IRAS 05491+0247 located in the L1617 cloud in Orion 
(D=460 pc). The blueshifted lobe is detected in the optical only 
from a distance of ~10" from the source. The first part of the 
beam is almost completely obscured by the molecular envelope, 
up to ~20", where it emerges from the cloud showing a chain of 
bright, equally spaced and well-collimated knots up to a distance 
of ~48" (knots E-L, following the nomenclature from Reipurth 
et al. 1997). These are followed by a series of more spaced and 
fainter knots up to 80" (M, N, O, P), and then the jet terminates 
with a strong bow-shock located at ~150" (knot V). The phys- 
ical properties of HH 111 have been investigated by several au- 
thors using both optical (e.g., Morse et al. 1993; Hartigan et al. 
1994) and near-infrared (e.g., Nisini et al. 2002) lines. Reipurth 
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Fig. 1. Position- Velocity diagrams of the forbidden line ratios used in the diagnostic technique for the HH 111 jet. From top to bottom panel: 
the [Sn]/16731 line intensity; the [S ii]/16731/6716 ratio increases with electron density, n e ; the [Nii]/[Oi] ratio mainly depends on the ionisation 
fraction, x e , and increases for increasing x e ; the [O i]/[S n] ratio increases for increasing temperature, TV 



et al. (1997) and Davis et al. (2001), instead, focused on the gas 
kinematics, through analysing [Sn] and H2 line profiles, and 
highlighted the presence of two velocity components peaking at 
— 75/-85 and ~-15/-30 kmr 1 . In a previous work (Podio et al. 
2006) a detailed analysis of the HH 111 physical structure was 
performed using several spectral tracers, which showed show- 
ing the density and temperature stratification in each spatially 
unresolved cooling region. Nevertheless, none of the previous 
studies investigated the variation of the physical parameters as a 
function of the gas velocity, i.e. in different velocity intervals. 

As we explain in Sect. 2 the data of HH 111 show partially 
velocity resolved line profiles. This is confirmed by the position- 
velocity (PV) diagrams of the computed line ratios in Fig. 1, 
which clearly show different value of the ratios in the two ve- 
locity intervals in the jet beam (knots from E to L). A compar- 
ison of the values of the ratios integrated over the two consid- 
ered velocity intervals (LVI and HVI) is also shown in Fig. 2. 
Since these ratios depend on the gas physical conditions, they 
give a qualitative idea of the variations of electron density, tem- 
perature, and ionisation fraction along the jet and with the gas 
velocity. The PV of the ratio between the sulphur doublet lines 
([S ii]/1673 1/6716), for example, shows that n f peaks at the posi- 
tions of the brightest knots and at high velocities. The [N n]/[0 1] 
and [Oi]/[Sii] ratios, on the other hand, indicate a peak of the 
excitation conditions (x e , T c ) in the high velocity interval of knot 
L. Interestingly, HST high-angular resolution images of the jet 
clearly show that HH 1 1 1 is a chain of shock working surfaces 
where knot L has the clearest "bow" morphology (Reipurth et al. 
1997) and is the knot with the largest [S 11] and H2 line profiles 
(Reipurth et al. 1997; Davis et al. 2001) and the largest shock 
velocity (Hartigan et al. 2001) in the jet beam. Our PV diagrams 
confirm this picture highlighting the shock velocity structure. 

By integrating the line fluxes over the seeing FWHM and 
over the LVI and HVI and then by applying the BE technique 
we obtain the variations of the physical parameters both along 
the jet and in the two velocity intervals for knots E-L (see Fig. 3). 



The derived parameters are in agreement with previous estimates 
(Podio et al. 2006): n e varies between 10 2 and ~2.2 10 3 cm -3 , 
x e goes from a few percent to ~0.3, T c ~0.7-1.3 10 4 K, and n# 
ranges between ~2 10 3 and ~3 10 4 cirT 3 . With respect to previ- 
ous analyses, however, the increased spatial sampling (more than 
twice with respect to previous studies) and the velocity informa- 
tion highlight several shocks features, such as n f maxima in the 
HVI of the brightest knots, higher values of T e in the HVI along 
the full jet length, and, finally, a strong and steep increase of x c 
and T e in the HVI of knot L. Note that very weak emission in the 
[Oiii],LI4959,5007 lines was detected by Morse et al. (1993) in 
this knot. The presence of oxygen in [O m] form is not accounted 
for in the adopted diagnostic technique. Nevertheless, since the 
emission in these lines is very weak, we tentatively applied the 
diagnostic to knot L and the obtained high values of x e and 
T e are in agreement with the detection of such high-excitation 
lines. These results are an important observational validation of 
the shock velocity structure predicted by models (Hartigan et al. 
1987). 

We do not apply the diagnostic to knots M, N, O, and P be- 
cause of the very weak emission in the [N 11] lines. Nevertheless, 
the electron density can be inferred by integrating the [S 11] emis- 
sion over the knot spatial profile to increase the S/N, which is 
low at the position of these faint knots. We obtain n e values of 
about 150 cm 3 , 200 cm 3 , 250 cm 3 , and 150 cm 3 in the HVI 
of knot M, N, O, and P and values of ~100 cirT 3 lower in the 
LVI for the same knots. The average values are in agreement 
with previous results from Podio et al. (2006). 

Previous works did not detect high excitation lines at the po- 
sition of the terminal bow shock (knot V) (Morse et al. 1993), 
thus we apply our diagnostic to this knot. Even if two veloc- 
ity component has been detected in the [S 11] lines by Reipurth 
et al. (1997) the PVs in Fig. 1 do not show different values of 
the ratios in the two velocity intervals. This can be due to the 
low spectral resolution combined with the much lower inten- 
sity of knot V with respect to the knots in the jet beam (knot 
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Fig. 2. Variation of the forbidden line ratios used in the diagnostic tech- 
nique for the HH 1 1 1 jet as a function of distance from the source and in 
two velocity intervals: the High Velocity Interval (HVI, red points/line, 
v<-100 kms~') and the Low Velocity Interval (LVI, black points/line, 
v>-100 km s~'). In contrast, in the terminal bow (knot V) the ratios has 
been computed by integrating the fluxes over the line spectral profile 
thus there is no information on the velocity intervals. From top to bot- 
tom panel: intensity profiles of the optical lines, [S n]/1673 1/67 16 ratios, 
[Nii]/[Oi] ratios, and [O i]/[S n] ratios. 



E-L). Thus we integrate the line fluxes over their spectral pro- 
files obtaining the trend shown in Fig. 3. The electron and total 
density are of ~400-550 cirr 3 and ~2-5 10 3 cm 4 , respectively. 
The ionisation fraction and the temperature show an unexpected 
anti-correlation and varies between ~0.23 and ~0.07, and ~0.75 
10 4 K and ~0.85 1 4 K moving from the shock front toward 
the source. Note that despite the high shock velocity of knot V 
(Hartigan et al. 2001) we did not infer high values of x e and T e 
as for knot L. 

4.2. The Sz68 jet 

The Sz68 jet has never been analysed through spectral diagnostic 
techniques. This jet is located in the B228 molecular cloud be- 
tween the two emission-line stars Sz68 and Sz69. The collimated 
chain of [S n] shock-excited blueshifted emission extending up 
to 34" from Sz68 at a position angle of 135° was detected by 
Heyer & Graham (1989). They report three knots located at 21", 
28", and 34" from the source, the brightest ones being the first 
knot, A, and the last one, B. Comparing the position of the knots 
in the spectra of Heyer & Graham (1989) with the positions in 



E F G H 11 J K 1 L V 




distance from source (arcsec) 



Fig. 3. Variation of the physical parameters for the HH 1 1 1 jet as a func- 
tion of distance from the source and in two velocity intervals: the High 
Velocity Interval (HVI, red points/line, v<-100 kms~') and the Low 
Velocity Interval (LVI, black points/line, v>-100 kms~'). In contrast, 
in the terminal bow (knot V) the parameters are derived integrating the 
fluxes over the line spectral profile thus there is no information on the 
velocity intervals. From top to bottom panel: intensity profiles of the 
optical lines, the electron density, n e , in units of 10 3 cm" 3 , the ionisa- 
tion fraction, x e , the temperature, T e , in units of 10 4 K, and the total 
density, n H , in units of 10 4 cm 4 . 



our EFOSC2 data (see upper panel of Fig. 4 for the knots spatial 
profile), we derive an estimate of the knot proper motions. This 
indicates a tangential velocity of ~600 km s for knot A, -460 
kmr 1 for the intermediate knot, and ~330 kms 1 for knot B, 
which combined with the radial velocity corrected for the source 
heliocentric velocity ( — 55 kms" 1 for knot A and — 29 kms 1 
for knot B, Heyer & Graham 1989) gives an upper limit for the 
inclination angle out of the plane of the sky of ~5°. 

We applied the BE technique to the line fluxes measured 
along the three knots by integrating all over the line spectral pro- 
files. Since this object is very faint, the error-bars are much larger 
than for the other jets. The inferred electron density is quite low 
varying between 200 ctrT 3 and 1.1 10 3 cirT 3 . These values are 
slightly lower with respect to those derived by Heyer & Graham 
(1989). The ionisation fraction is 0.05-0.25 and the temperature 
ranges between 0.6-1.1 10 4 K, with a maximum at the position 
of the intermediate knot, where there is a peak of n e , as well. 
Finally the hydrogen total density is 2-6 10 3 cm 4 . 

The determination of the jet physical conditions, as well as 
the jet tangential velocity and inclination angle allows us to esti- 
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distance from source (arcsec) 

Fig. 4. Variation of the physical parameters for the Sz68 jet as a function 
of the distance from the source. From top to bottom panel: intensity 
profiles of the optical lines, the electron density, n e , in units of 10 3 cm" 3 , 
the ionisation fraction, x c , the temperature, T c , in units of 10 4 K, and the 
total density, n H , in units of 10 4 cirr 3 . 



mate the mass loss rate as Mj et = \irnn riHTtrj vj, where //=1.41 
is the mean atomic weight, m# the proton mass, tin the hydrogen 
density and r j and v j, respectively, the jet radius and velocity. By 
assuming a typical jet radius of ~ 100-300 AU we obtain a mass 
loss rate between a few 10~ 8 and a few 10~ 7 M©yr~'. These val- 
ues are of the same order of magnitude of Mj et found for typical 
HH jets (Hartigan et al. 1994; Bacciotti & Eisloffel 1999; Podio 
et al. 2006). 



5. Dust Reprocessing in the Jets 

Refractory species, such as Ca, Fe, Ni, Cr, Ti, Mg, Si are of- 
ten depleted in the interstellar medium because their atoms are 
locked onto dust grains (Savage & Sembach 1996). The calcium 
gas-phase abundance, for example, is estimated to be ~2 10~ 8 in 
Orion (Baldwin et al. 1991), i.e. two orders of magnitude lower 
than its Solar abundance (-2.04 10~ 6 , Asplund et al. 2005), and 
up to 4 order of magnitude lower than Solar in f Oph (Savage 
& Sembach 1996). On the other hand, the passage of shocks can 
partially or completely destroy the dust releasing the refractory 
atoms into the gas cloud. 

Theoretical models of dust reprocessing in shocks, such as 
those driven by supernova explosions or those occurring in stel- 
lar jets, estimated the dust destruction rate by considering the 
different processes at work, i.e. thermal and inertial sputtering, 



photoevaporation and shattering. They showed that the relative 
importance of these processes, and thus the efficiency of the 
shock in reprocessing the dust, strongly depend on a number of 
parameters such as the shock velocity, the gas pre-shock density, 
the intensity of the magnetic field, the size and structure of the 
dust grains (e.g., Jones et al. 1994; Jones 2000; Draine 2003; 
Guillet et al. 2009). 

From an observational point of view an estimate of the de- 
pletion of refractory species with respect to their Solar abun- 
dance can be used to gauge the amount of dust grains in stellar 
jet beams and thus to evaluate the efficiency of these mild shocks 
in destroying the dust. This in turn can help constraining the lo- 
cation and size of the region around the star from where the jet 
originates. 

In Sect. 5. 1 we explain in detail the method used to estimate 
the abundance of calcium in the gas-phase, while in Sect. 5.2 we 
present the results obtained and in Sect. 5.3 we compare them 
with the predictions of theoretical models. 

5.1. Observational estimates of calcium gas-phase 
abundance in stellar jets: the method 

To estimate the depletion of a given refractory specie, in our 
case calcium, with respect to its Solar abundance, we com- 
pare "expected" and observed ratios between refractory and non- 
refractory species known to be present in Solar abundance in the 
ISM (e.g. Savage & Sembach 1996). 

The "expected" ratios can be computed through two different 
methods: 

(i) by using shock models which predict the gas physical 
conditions and thus the emission in different lines behind the 
shock front. This method was used by Mouri & Taniguchi 
(2000) and, in part, by Nisini et al. (2002) to estimate Fe 
gas-phase abundance in HH jets. They found that Fe is 
strongly depleted in these objects, up to 80% with respect to 
its Solar abundance; 

(ii) by using the gas physical conditions inferred through spec- 
tral diagnostic tools. In this case the estimate does not de- 
pend upon assumptions of the shock characteristics (veloc- 
ity, pre-shock density, magnetic field intensity) but depends 
on the accuracy of the determinations of the gas parameters. 
Beck-Winchatz et al. (1994, 1996) used this method to esti- 
mate Fe and Ni gas-phase abundance in HH jets. They found 
that Fe is not, or only partially, depleted in HH jets, while 
Ni is surprisingly over-abundant. The method was then used 
by Nisini et al. (2002, 2005) and Podio et al. (2006) to esti- 
mate the gas-phase abundance of both Fe and Ca obtaining 
a strong depletion (up to 90%) of these species which, how- 
ever, may vary along the jet and decreases down to 0% in 
some knots. 

The main assumption of method (ii) is that the considered 
lines from refractory and non-refractory species are emitted in 
the same region behind the shock front, in such a way that the 
filling factor and the gas physical conditions are the same. This 
has been verified by Nisini et al. (2005) (see their Fig. 3 and 
corresponding discussion). 

Following Nisini et al. (2005) and Podio et al. (2006), here 
we use the [Can]/l7291/[Sii]/16731 ratio to estimate the Ca gas- 
phase abundance in our sample of HH jets. These lines have 
similar excitation temperatures (~2 10 4 K for [Can]/I7291 and 
~2.210 4 K for [Sn]/I6731) and, despite the critical density 
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of the calcium line being larger than that of the sulphur line 
(n cr ([Sn])~2.5 10 4 cm- 3 , n cr ([Can])~5 10 7 cm" 3 ), at the typi- 
cal low density of HH jets the emission of the two lines peaks 
at the same distance from the shock front. In the same region 
the [S n], [Oi], and [Nn] lines which are used to derive the gas 
physical conditions along the jet are also excited. Thus the esti- 
mated parameters can be used consistently to model the theoret- 
ical [Can]/[S n] ratio. 

Following Hartigan et al. (2004) we assume that the [Ca n] 
lines are collisionally excited and compute the level populations. 
As in previous works (Nisini et al. 2005; Podio et al. 2006) we 
assume that there is no calcium in the form of Ca°, because its 
ionization potential is very low, ~6.1 eV. On the other hand, the 
ionization potential of Can is also quite low, ~11.9eV, i.e. lower 
than that for hydrogen (~13.6eV). In Nisini et al. (2005) and 
Podio et al. (2006) we computed the [Can]/[Sii] ratios by as- 
suming that all calcium is in the form of Ca + . This is a reason- 
able approximation in the case of poorly ionized objects such as 
HH 1 and HH 111 (x e <0.3). However, the presence of Ca in the 
form of Ca ++ cannot be neglected when analysing higher exci- 
tation jets such as HH 83 and HH 24. The coefficients for col- 
lisional ionization and radiative and dielectronic recombination 
of Ca + and H° are very similar for temperatures between 8 10 3 
K and 3 10 4 K (at T~2 10 4 K, the recombination rates for Ca + 
andH are: a(Ca + )~3.9 10~ 13 cm 3 s _1 , or(H°)~2.5 10~ 13 cm 3 s" 1 ; 
while the coefficients for collisional ionization are: C(Ca + )~1.2 
10" 11 cm 3 s -1 , C(H°)~1.6 lO" 11 cm 3 s" 1 ). Therefore, considered 
that at these temperatures all calcium is in the form of Ca + and 
Ca ++ , in conditions of ionization equilibrium the ionization frac- 
tion of Ca + (defined as x(Ca + )=Ca ++ /Ca) and H° are almost the 
same: x(Ca + )~x e . In stellar jets, however, the ionization fraction 
of hydrogen is never at equilibrium with the local temperature 
because, as explained in Bacciotti & Eisloffel (1999), the gas 
is moving along the jet and the recombination timescale at the 
observed electron densities is of ~10 3 yr, i.e. of the same or- 
der of the jet dynamical time. Because of the similarity of the 
recombination and collisional ionization coefficients, however, 
Ca + and H° atoms can be thought to be undergoing the same 
processes when moving along the jet. Thus we assume that, at 
every position along the jet, the Ca + ionization fraction is equal 
to the hydrogen ionization fraction, x e . Taking Ca and S gas- 
phase abundances equal to Solar ones (Asplund et al. 2005), we 
can finally compute the [Ca n]/[S n] theoretical ratios. 

Note that, in contrast with previous works (Nisini et al. 2005; 
Podio et al. 2006), here we use the [Can]^7291/[S ii]/16731 
ratio, instead of the [Can]/l(7291+7324)/[Sii]^(6716+6731), 
since the [Can]/!7324 line may be blended with the 
[On],U7318.6,7319.4 and the [On],U7329.9,7330.7 lines. As 
explained in Hartigan et al. (2004) the calcium line ratio 
[Ca n]/1729 1//17324 is expected to be ~ 1 .5 both in the case of col- 
lisional excitation and for fluorescent pumping from the ground 
state. Thus by measuring the ratio between the calcium lines we 
can check for the presence, if any, of [On] emission blended 
with the [Can]/17324 line. We find that this is not the case for 
the HH 1 and the HH 34 jets (Nisini et al. 2005; Podio et al. 
2006), but may be the case in the HH 1 1 1 , HH 83 and HH 24 jets 
which show higher excitation conditions in some of the knots 
(see Figs. 3, A. 3, and A.4). These jets, in fact, show a wider 
[Can]/17324 line in some knots indicating blending with the 
[On] lines. Moreover, in our EFOSC2 spectra, the [Can]/i7324 
line is at the border of the detector and hence the line flux can be 
affected by distortion and/or higher S/N. Thus we decided to use 
the [Caii]/l7291/[Sn]/16731 for all the jets in our sample. 



Table 1. Gas-phase abundance of calcium with respect to Solar 
Abundances. 



Jet 


[Ca] ;nj /[Ca] 5 ' ; fl / 


HH 111 


0.21-0.52 


HH 1 


0.25-0.79 


HH 83 


0.53-1.57* 


HH 24 A 


0.40-0.69* 


HH 24 E 


0.62-0.85* 


HH 24 C 


0.52-0.82* 


HH 34* 


0.3-1 


Orion Cloud c 


0.01 



"(Ca/H) Solar abundance from Asplund et al. (2005): 2.04 10" 6 
*from Podio et al. (2006) 

c (Ca/H) abundance estimated in the interstellar medium in Orion by 
Baldwin etal. (1991): 2 10- 8 

*the values of [Ca.] gas /[Ca.] s „i ar for HH 83 and HH 24 have been esti- 
mated assuming A v =0, thus they are upper limits to the real values. 



For the HH 111 and the HH 1 jet we corrected the observed 
[Can]/[Sn] ratios for reddening following the procedure ex- 
plained in Sect. 3. For the HH 83 and the HH 24 jet we have 
no estimates of the visual extinction, thus we could not correct 
the Ca + and S + line fluxes for reddening and the observed ra- 
tios (computed assuming A v =0) are actually upper limits to the 
real values (the latter decrease for increasing Ay). This means 
that the calcium depletion inferred from Figs. 7 and 8 is a lower 
limit. 

The errors affecting the observed ratios are due to the signal- 
to-noise over the measured line fluxes. For HH 1 and HH 111 
we consider also the errors due to the uncertainty in estimating 
Ay. 

The errors on the theoretical ratios are due to the uncertainty 
of the estimated physical parameters and are computed by vary- 
ing n e , x e , and T e between their minimum and maximum values. 
The influence of these errors is represented in the form of a green 
stripe in Figs. 5, 6, 7, and 8. Since, as explained in Sect. 3, the 
choice of the abundance set may affect the absolute values of x e 
and T e , we also computed the variations of [Can]/[Sii] when 
adopting the alternative abundance set estimated by Esteban et 
al. (2004). The difference between the [Can]/[Sn] ratios ob- 
tained with the two sets is illustrated by the width of the beige 
stripes in Figs. 5, 6, 7, and 8. 

5.2. Results: Ca depletion along HH jets. 

The comparison between the observed and the predicted 
[Ca n]/[S n] ratios is shown in Figs. 5, 6, 7, and 8 and the inferred 
values of the calcium gas-phase abundance are summarized in 
Table 1. Note that [Ca] gas /[Ca] So i ar values in Table 1 and in the 
bottom panel of Fig. 5 are derived by using the [Can]/[S n] ratios 
predicted assuming Solar elemental abundances from Asplund et 
al. (2005) (illustrated by a dotted line in the figures). 

Along the HH 111 and HH 1 jets the observed 
[Ca n]T729 1/[S n]i673 1 ratios are lower than the predicted ones, 
i.e. calcium is depleted with respect to its Solar abundance, the 
depletion being between 20% and 80%. For these jets the results 
agree with what has been found in Nisini et al. (2005) and Podio 
et al. (2006), the slight difference being explained by the fact that 
in those works the presence of calcium in the form of Ca ++ was 
neglected and a different set of abundances was assumed. 

On the contrary, along HH 83 the dust grains appear to have 
been destroyed substantially, as the derived Ca abundance is 
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Fig. 5. Comparison between observed (diamonds) and predicted 
(dashed lines) [Can]A7291/[S n]i6731 along the HH 1 1 1 jet and in two 
velocity intervals: the HVI (v<-100 km s -1 , top panel) and the LVI (v>- 
100 km s , bottom panel). The dotted line refers to the predicted ratios 
computed assuming Solar abundances from Asplund et al. (2005) and 
the green stripe shows the variation of the ratios due to the errors on the 
inferred physical parameters. The beige stripe, on the contrary, shows 
the variations of the predicted ratios when assuming non-solar abun- 
dances from Esteban et al. (2004). In any case, the difference between 
observed and predicted ratios indicates a depletion of the Ca gas-phase 
abundance, which is quantified in the bottom panel using the predicted 
ratios computed assuming Solar abundances (dotted lines). 



equal or even larger than Solar in the first knots and is depleted 
up to 47% at distances larger than 25" from the source. In the 
case of HH 83, however, we did not correct the fluxes for ex- 
tinction, thus the obtained values of Ca gas-phase abundance are 
upper limits. It can be shown that with an Ay=5 the Ca gas-phase 
abundance is between 0.36 and 1. In any case the low depletion 
in this jet appear to be consistent with the jet high excitation. 

Finally, in HH 24 we derive the [Can]^7291/[Sn]^6731 ra- 
tio only in a few of the observed knots because for the fainter 
ones the Ca + emission is below our assumed S/N threshold 
(equal to three times the background noise). Fig. 8 shows that 
the depletion along the jet, show no evident trend: calcium is 
depleted by 48%-18% in HH 24 C, by ~38%-15% in HH 24 E 
and by ~60%-31% in HH 24 A. The lack of a trend along the 
jet reflects the situation found for the physical parameters (see 
Fig. A. 4) and supports the idea that HH 24 C, E, and A may 
not belong to a single HH jet (Eisloffel & Mundt 1997). Note 
that also in the case of HH 24 we did not correct the observed 
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Fig. 6. Comparison between observed (diamonds) and predicted 
(dashed line) [Can]/17291/[S n]i6731 along the HH 1 jet. The green 
stripe and the beige stripe show the variation of the predicted ratios, 
as explained in Fig. 5. The difference between observed and predicted 
ratios indicates a depletion of the Ca gas-phase abundance. 




15 20 25 30 

distance from source (arcsec) 

Fig. 7. Comparison between observed (diamonds) and predicted 
(dashed line) [Can]/I7291/[S HR6731 along the HH 83 jet. The ob- 
served ratios are not corrected for reddening (Ay=0), and are thus upper 
limits. The green stripe and the beige stripe show the variation of the 
predicted ratios, as explained in Fig. 5. 



[Ca n]/[S n] ratios for extinction thus the inferred Ca depletion 
is a lower limit. 

The width of the beige stripes in Figs. 5, 6, 7, and 8 show 
the variation of the theoretical [Ca n]/[S n] ratios when adopting 
abundances from Esteban et al. (2004). Note that also in this case 
calcium turns out to be largely depleted in most of the knots, 
although the inferred Ca depletion is lower. 
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Fig. 8. Comparison between observed (diamonds) and predicted 
(dashed line) [Can]^7291/[S n]^6731 along the HH 24 jet. The ob- 
served ratios are not corrected for reddening (Ay=0) and thus they are 
an upper limits. The green stripe and the beige stripe show the variation 
of the predicted ratios, as explained in Fig. 5. 



5.3. Discussion. 

In general, we find an amount of Ca depletion comparable to 
the Fe depletion estimated in other works (Beck-Winchatz et al. 
1996; Bohm & Matt 2001; Mouri & Taniguchi 2000; Nisini et 
al. 2002, 2005; Podio et al. 2006; Garcia Lopez et al. in prepara- 
tion). Together the results obtained for the two different refrac- 
tory species indicate that, despite the presence of shocks, there 
are still dust grains in the jets. However, partial or total destruc- 
tion of dust grains can be inferred by comparing the level of 
calcium depletion in the jets with that of the diffuse interstellar 
medium in Orion (see Tab. 1). 

According to the most recent models of dust reprocessing 
from Guillet et al. (2009) the dust destruction rate in J-shocks 
with velocities lower than 50 km s~' and pre-shock densities 
from 10 4 to 10 6 cirT 3 is around a few percent. Previous models 
by Jones et al. (1994) predict a destruction rate for silicate grains 
as high as 50% for shock velocities of ~ 100 km s" 1 and low val- 

— 1/2 

ues of the b parameter (b — So« , where Z?o is the intensity 
of the magnetic field and «o is the pre-shock density). Finally, 
up to 90% of Fe and 60%-70% of Mg and Si can be released 
in gaseous form in C-shocks with velocities of only ~45 km s" 1 
(May et al. 2000). It is important to stress that none of the above 
mentioned models has been tested with grains containing atoms 
of Ca and that these predictions are highly dependent on dust 
size, structure, and composition. Different results, for example, 
are retrieved when considering silicate or graphite grains, and if 
the grains are solid or porous (Jones et al. 1994; Guillet et al. 
2009). 

In any case, none of the present models predict a total dust 
grain destruction in the mild shocks that occurr in the jet beam 
where the material interacts with gas already put into motion by 
the passage of previous fronts (v s <100 km s" 1 in the jet beam, 
Hartigan et al. 2001). This discrepancy between the predictions 
of shock models and the measured Ca gas-phase abundance is 
easily solved if we assume that the dust has been reprocessed in 



the disk or circumstellar region before being extracted and ac- 
celerated into jets. According to recent theoretical models and 
observations, the accretion disk in Herbig and Classical T Tauri 
Stars is populated by dust grains only beyond the so-called dust 
evaporation radius, R evp (e.g., Isella & Natta 2005). One can 
thus infer that the material in the jet which is coming from a re- 
gion smaller than R evp would not contain grains, as these have 
been destroyed by the stellar radiation. The dust evaporation ra- 
dius is located at ~0. 1-1 .5 AU in Herbig stars (Eisner et al. 2007) 
and at -0.05-0.3 AU in T Tauri stars (Akeson et al. 2005), while 
there are no estimates in younger and more embedded sources 
such as the emitting sources of the HH jets in our sample. 

On the other hand, the depletion is observed at various lev- 
els in the beam, and to explain the presence of dust grains we 
have to assume that part of the material in the jet is coming from 
a region of the disk which extends beyond R evp . Note that we 
are not observing ambient dust since this has been swept-out or 
destroyed by the passage of the leading bow-shock, which has 
larger shock velocities (v s up to ~150 km s" 1 in the bow of HH 
1, Bally et al. 2002) and the dust reformation time-scale is much 
larger than the dynamical timescale of the jet, which is of 10 3 - 
10 4 yr. Moreover, the rise in gas-phase calcium with distance 
from the source along the HH 111 and HH 1 jets argues against 
major entrainment of matter along the jet, as would be suggested 
by models of prompt and/or turbulent entrainment (Masson & 
Chernin 1993; Stahler 1994). 

Interestingly, in HH 111, for which we could perform the 
analysis in the two velocity intervals, the Ca depletion is at a 
minimum in the HVI of the bow-shaped knot L (see HST im- 
ages from Reipurth et al. 1997) which is also the one with the 
higher shock velocity (v,~60 kms~', Hartigan et al. 2001). This 
result suggests that the efficiency of shocks in destroying dust 
grains strongly depends on the shock velocity, in agreement with 
the predictions of theoretical models (Jones et al. 1994; Draine 
2003; Guillet et al. 2009). Our estimate also indicates that this 
shock should be able to destroy at least 40% of the dust grains 
present at this position. 

The proposed proportionality of destruction efficiency and 
shock velocity is also supported by our finding that calcium is 
less or not depleted along HH 83 and HH 24. These jets, in fact, 
show higher excitation conditions in comparison to HH 1 and 
HH 111 suggesting that their gas has undergone stronger internal 
shocks. 

6. Summary and Conclusions 

In this paper we analyse partially velocity resolved optical spec- 
tra of a sample of HH jets located in the Orion and in the Lupus 
Molecular Cloud (HH 111, HH 1/2, HH 83, HH 24 M/A/E/C, 
and Sz68), to investigate primarily the content of Ca in gas-phase 
and, in turn, the dust reprocessing in jets. To this aim we first ap- 
ply spectral diagnostics (Bacciotti & Eisloffel 1999) to derive 
the physical conditions of the gas in the jets from the line ratios 
of the detected forbidden lines, i.e. the values of the electron and 
total density, the ionisation fraction, and the temperature. These 
parameters are then used to estimate the expected Ca abundance. 

We find that in the studied jets the electron density ranges be- 
tween 50 and 4 10 3 cm" 3 . The ionisation fraction varies between 
a few percent in low excited knots where the [N n] emission is 
very faint and -0.4-0.7 in high excitation knots with strong [N n] 
emission. The temperature follows the same trend as the ioni- 
sation fraction, being ~10 4 K in the weakly ionised knots and 
~1.5-3 10 4 K in the highly ionised ones, in agreement with pre- 
dictions from shock models (e.g., Hartigan et al. 1987). Finally, 
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the total density derived by combining the retrieved values of n e 
and x e varies between ~100 and 6 10 4 cm 4 . 

The most interesting results are obtained for the HH 1 1 1 jet 
for which we derive the values of the physical parameters sepa- 
rately in two velocity intervals, Low (v>-100 kms 1 ) and High 
(v<-100 kms 1 ) (LVI and HVI). This analysis highlights many 
shock signatures such as n e maxima in the HVI of the brightest 
knots and a sharp increase of x e and T e in the HVI of the strong 
knot L, which has a clear bow morphology in the HST images 
of Reipurth et al. (1997) and is the knot with the highest shock 
velocity in the jet beam (Hartigan et al. 2001). 

We also analysed for the first time the faint jet powered by 
the young source Sz68. In this case we obtain low values of the 
electron and total density (n c ~200-10 3 ctrT 3 and n H ~2-6 10 3 
cm -3 ) and intermediate values of the ionisation fraction and the 
temperature (x e ~0.05-0.25 and T e ~0.8 10 4 K). For this jet we 
also measured proper motions of ~300-600 km s _1 from which, 
assuming a jet radius of ~100-300 AUs we derived a mass loss 
rate ranging between a few 10~ 8 and a few 1O~ 7 M yr _1 , i.e. of 
the same order of magnitude of Mj et found for typical HH jets 
(Hartigan et al. 1994; Bacciotti & Eisloffel 1999; Podio et al. 
2006). 

Using the inferred physical parameters we derive an es- 
timate for the dust content along the jets by comparing ob- 
served and predicted ratios between calcium and sulphur lines. 
The derived Ca gas-phase abundance is lower than Solar 
([Ca] gas /[Ca]soiar~0.2-0.S5) in all the jets, with the exception of 
HH 83 to which, however, no dereddening has been applied, thus 
yielding upper limits to the Ca gas-phase abundance. In general, 
the measured Ca depletion is lower with respect to the inter- 
stellar medium ([Ca] 4aj /[Ca] lSo ; ar ~0.01 in Orion, Baldwin et al. 
1991). This result suggests that mild shocks occurring along the 
jets are partially destroying dust grains, thus releasing calcium 
atoms into gaseous form. Interestingly, a higher degree of deple- 
tion is estimated in the low-excited jets (i.e. HH 111 and HH 1), 
while small or no depletion is found along HH 83 and HH 24 
which show higher ionization fractions and temperatures. This 
seems to indicate that the efficiency of shocks in reprocessing 
dust grains strongly increases with shock velocity, as expected 
from theoretical models (e.g., Jones et al. 1994; Guillet et al. 
2009). This idea is also supported by the finding that along HH 
111 the depletion is minimum in the HVI of knot L, which has 
the largest shock velocity (Hartigan et al. 2001). 

The measured values of Ca gas-phase abundance, however, 
are too large to be explained by dust reprocessing in shocks 
which are slower than 100 kms -1 (Jones et al. 1994; May et 
al. 2000; Guillet et al. 2009). This discrepancy between observa- 
tional results and predictions of theoretical models may indicate 
that the material in the jet is extracted from a wide region of the 
disk that includes both the inner region close to the star, inside 
the so-called dust evaporation radius, R evp , where the dust is de- 
stroyed by the stellar radiation (Isella & Natta 2005), and the 
region beyond R evp from where the dust is lifted and accelerated 
in the flow. Modeling, however, is required to test this idea. 
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Fig. A.l. Variation of the physical parameters for the HH 1 jet beam as 
a function of the distance from the source. From top to bottom panel: 
intensity profiles of the optical lines, the electron density, n f , in units of 
10 3 cm" 3 , the ionisation fraction, x c , the temperature, T„, in units of 10 4 
K, and the total density, n H , in units of 10 4 cm" 3 . 

Appendix A: Physical structure of previously 
analysed jets 

In this appendix we present the results obtained by applying 
spectral diagnostics to HH 1 (jet and bow), HH 2, HH 83, and 
HH 24. These jets have been already analysed in previous papers 
with the same diagnostics, but are re-analysed here with higher 
spatial sampling (Ax~1.2") and an improved version of the di- 
agnostic code as explained in Sect. 3. The values presented are 
used for the investigation of the dust reprocessing in these jets 
(see Sect. 5). 

A. 1 . HH 1/2: the jet and its terminal bows 

The physical structure of the HH 1 jets and its terminal bows 
has been previously analysed by Solf & Bohm (1991), Nisini et 
al. (2005), and Bohm & Solf (1985), Solf et al. (1988), Eisloffel 
et al. (1994), Bally et al. (2002). In our observations the slit is 
aligned along the jet and it covers the fainter western/eastern part 
of HH 1 and HH 2 respectively, i.e. the knots B and G, which are 
located to the west of the HH 1 bow apex (knot F), and the knots 
L, J, G, B, T, and Q located to the east of the bright H and A 
knots in HH 2. 

The results obtained by applying the BE technique to the HH 
1 jet are shown in Fig. A.l. The values for the physical parame- 



ters are in a good agreement with the ones obtained in Nisini et 
al. (2005) but here the spatial sampling is two times better. 

We could not apply the diagnostic to the terminal bows. The 
BE technique, in fact, relies on the assumption of low excita- 
tion conditions, i.e. S is all ionised but there is no S ++ , and oxy- 
gen and nitrogen are ionised at most once (Bacciotti & Eisloffel 
1999). This assumption is satisfied when studying the jet beam 
where the gas interacts with material that was already put into 
motion from the passage of previous shocks, but may not be cor- 
rect when dealing with the terminal bows. This picture is con- 
firmed by proper motions studies by Bally et al. (2002) which 
indicate shock velocities lower than 30 km s" 1 in the jet and ve- 
locity variations of up to 100-200 km s" 1 in the bows. Moreover, 
Bohm & Solf (1985) and Solf et al. (1988) detected many high 
excitation lines in HH 1 and HH 2 such as ++ , S ++ , and ArlV. 

Even if in our case the slit is covering the lateral part of the 
bows the excitation conditions may still be too high. This is indi- 
cated by the detection of the high excitation [Arm]/!7135.8 line 
in our spectra and by the fact that the [N n] lines, which are faint 
in the jet beam, show strong emission in the bows, comparable 
to the [S n] lines (see the upper panel of Fig. A. 2). 

In order to compare the physical conditions in the jet and 
in the bows we computed the electron density, n c , and the 
[Nii]/[Oi] ratio along all the slit length. The electron den- 
sity does not show large differences. It varies between 0.05- 
4 10 3 cm" 3 along all the jet. While in the jet beam n e is de- 
creasing with distance from the source, both in the HH 1 and in 
HH 2 bows, on the contrary, there is a decreasing trend going 
from the shock apex toward the source. This is expected behind 
a shock front and confirms the results found by Bohm & Solf 
(1985) and Solf et al. (1988). Our values of n e are lower than 
those inferred by Bohm & Solf (1985), however, because of the 
different slit alignment (on the brightest spots at the apex of HH 
1 and HH 2 for Bohm & Solf 1985, and on the bows wings in 
our case) and show that the density in the two bows is maximum 
at the apex and is fading towards the wings. The [Nii]/[Oi] ratio 
is a good indicator of the excitation conditions and, in particular, 
of the ionisation state of the gas. Fig. A. 2 shows that [Nii]/[Oi] 
is <1 in the jet beam and in knot L, while it is >1 in the HH 1 
and HH 2 bows, indicating that the excitation is much higher in 
the bows, and a higher value of x e is expected. 

A.2. HH 83: the physical structure of the inner knots 

The physical structure of the HH 83 jet has been already derived 
in Podio et al. (2006). Thanks to the high S/N of these data, how- 
ever, we obtained a sampling which is four times larger with re- 
spect to previous results. Moreover the good quality of the data, 
which allowed us to properly subtract the continuum emission 
from the reflection nebula Re 17 (Rolph et al. 1990), and the 
use of the "improved" diagnostic code allowed us to estimate 
the gas physical conditions in the inner part of the jet, where the 
emission in the [N n] lines is comparable to the [S n] and [O i] 
emission. 

The derived parameters indicate that the excitation condi- 
tions are very high in the knots near to the source (knots from 
A to D) with n e ~400-700 cm" 3 , and high values of the ionisa- 
tion fraction and temperature (x e ~0.4-0.7, and T e ~1.5-2 10 4 K). 

A.3. The HH 24 jets 

The physical parameters along the jets HH 24 C, E, and A, has 
been already derived in Bacciotti & Eisloffel (1999) and Podio 



Linda Podio et al.: Physical Structure and Dust Reprocessing in a sample of HH Jets, Online Material p 2 



HH 2 BOW 
B1 

B2 T Q 



3 — 



2 — 



m 



Ft 



1.5 — 

1 .0 ; -J 

0.5 -w- 



HH 1 JET 
G FE B A 



HH 1 BOW 
G 




-100 



-80 



-60 -40 10 30 

distance from source (arcsec) 



50 



70 



Fig. A.2. Analysis of the excitation conditions in the HH 1 jet and its terminal bows HH 1 and HH 2 as a function of the distance from the source. 
From top to bottom panel: intensity profiles of the optical lines, the electron density, n e , in units of 10 3 cm 4 , and the [Nii]/[Oi] ratio. 




Fig. A.4. Variation of the physical parameters for the HH 24 C/E jet as a function of the distance from the source. From top to bottom panel: 
intensity profiles of the optical lines, the electron density, n e , in units of 10 3 cm" 3 , the ionisation fraction, x e , the temperature, T e , in units of 10 4 
K, and the total density, n H , in units of 10 4 cm 4 . 
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Fig. A.3. Variation of the physical parameters for the HH 83 jet as a 
function of the distance from the source. From top to bottom panel: 
intensity profiles of the optical lines, the electron density, n e , in units of 
10 3 cm" 3 , the ionisation fraction, x c , the temperature, T„, in units of 10 4 
K, and the total density, n H , in units of 10 3 ever 3 . 

et al. (2006). In our observations the slit has been aligned along 
the axis HH 24 M-A-E and thus only partially covers the knots 
of the HH 24 C jet up to 40" (knot C6). This is why the line 
profiles in the upper panel of Fig. A.4 show fainter emission in 
the knots of group C, contrary to what was found in Podio et al. 
(2006), where the slit was aligned along the HH 24 C jet. The 
variation of the physical parameters obtained by applying the 
BE technique is shown in Fig. A.4. The sampling is improved 
of around one third with respect to previous analyses (Bacciotti 
& Eisloffel 1999; Podio et al. 2006) allowing us to highlight the 
different excitation conditions in the various groups of knots de- 
tected in the HH 24 complex (HH 24 A, HH 24 M/E, and HH 
24 C) and supporting the idea that these knots may belong to 
different jets. 



